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What are catamorphisms? 

• A family of functions: often called fold operators. 

• They replace data constructors with arbitrary functions. 

Catamorphisms are useful: 

• they are elegant means of expressing algorithms; 

• they facilitate reasoning about programs; 

• they support program optmization (e.g., deforestation). 
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The list catamorphism 

Similar to ML’s fold and Haskell’s foldr. 
datatype List(a) = Nil | Cons of a x List(a) 

fun cata(n,c) Nil = n 

I cata(n,c) (Cons(a,r)) = c( a, cata(n,c) r) 

Function c replaces Cons, value n replaces Nil: 
cata(n,c) (Cons(1,Cons(2,Cons(3,Nil)))) = c(1 ,c(2,c(3,n))) 
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Example 




fun filter(p) x = 

cata( Nil, A,(a,r). if p(a) then Cons(a,r) else r) x 


filter(even) (Cons(1 ,Cons(2,Cons(3,Cons(4,Nil))))) 
= c(1 ,c(2,c(3,c(4,n)))) 

= Cons(2,Cons(4,Nil)) 
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Catamorphisms can be easily defined for most first-order 
algebraic datatype. 

datatype Term = Const of int 
Var of string 
Appl of Term x Term 
Abs of string x Term 

(kx. x) 5 -► Appl(Abs(“x”,Var “x”),Const 5) 

fun cataT(c,v,p,a) (Const n) = c n 
I cataT(c,v,p,a) (Var s) =vs 
I cataT(c,v,p,a) (Appl(x,y)) = p( cataT(c,v,p,a) x, 

cataT(c,v,p,a) y) 

I cataT(c,v,p,a) (Abs(s,e)) = a( s, cataT(c,v,p,a) e ) 

v_ J 


3 & Sheard 




-6- 









( \ 

Higher-order datatypes 

No problem, unless they contain a contravariant recursion : 

datatype T = C of int | D of T -»T 

Why bother? There are 3 examples: 

• higher-order abstract syntax; 

• circular lists; 

• graphs. 
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Higher-order abstract syntax 

Pfenning & Elliot [PLDI’88], Wand [POPL’93], Despeyroux & Hirschowitz [ LPAR’94\ 

Closed terms of ?i-calculus: 

datatype Term = Const of int 

I Appl of Term x Term 
I Abs of Term -> Term 

(Xx. x x) 5 -► Appl(Abs(fn x => Appl(x,x)),Const 5) 
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Evaluator over higher-order syntax 

datatype Term = Const of int 

I Appl of Term x Term 
I Abs of Term Term 

datatype Value = Num of int | Fun of Value->Value 

eval: Term Value 

fun eval(Const n) = Num n 
I eval(Appl(f,e)) = let val Fun(g) = eval(f) 
in g(eval(e)) end 
I eval(Abs f) = Fun( ...) 

eval( Appl(Abs(fn x => x),Const 1)) = Num 1 


Circular Lists 

datatype Clist(a) = Nil 

I Cons of ax Clist(a) 

I Rec of Clist(a) Clist(a) 

fun head(Cons(a,r)) = a 
| head(Recf) = head(f(Rec f)) 

The circular list [1,2,1,2,1,2,...] is: a ImE? 

Rec( fn x => Cons(1 ,Cons(2,x))) 

[1,2,3,4,...] is: 

Rec( fn x => Cons( 1, map(fn z => z+1) x )) 
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Graphs 

datatype Graph(a) = Node of ax List(Graph(a)) 

| Rec of Graph(a) -> Graph(a) 

| Share of (Graph(a) Graph(a)) x Graph(a) 


Share( fn x => T(x), e ) 

= let val x = e 
in T(x) end 

fun root(Node(a,r)) = a 
| root(Recf) = root(f(Rec f)) 

| root(Share(f,e)) = root(f e) 

Rec(fn x => Share(fn z => Node(0,[z,Rec(fn y => Node(1 ,[y,z]))]), 
Node(2,[x]))) 
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An example of our technique 

• the recursive Term evaluator in the style of Paterson- 
Meij er-Hutton; 

• the problems of this method; 

• the recursive Term evaluator in our style; 

• the catamorphism over Term. 
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The X term evaluator (a la Paterson-Meijer-Hutton) 

datatype Term = Const of int 

I Appl of Term x Term 
I Abs of Term Term 

datatype Value = Num of int | Fun of Value-Walue 


fun eval(Const n) = Num n 
I eval(Appl(f,e))= let val Fun(g) = eval(f) 
in g(eval(e)) end 
I eval(Absf) = Fun( G ) 

Need to define G : Value -» Value 
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i 

eval(Abs f) = Fun( G ) reify 
Value 

fun eval(Const n) = Num n 
I eval(Appl(f,e))= let val Fun(g) = eval(f) in g(eval(e)) end 
I eval(Abs f) = Fun( eval ° f ° reify ) 
and reify(Num n) = Const n 
! reify(Fun g) = Abs( reify ° g ° eval) 
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Evaluating (Xx. x) 1 using the eval-reify approach 

eval( Appl(Abs(fn x => x),Const 1)) 

= let val Fun(g) = eval(Abs(fn x => x)) 
in g(eval(Const 1)) end 

= let val Fun(g) = Fun(fn x => eval(reify(x))) 
in g(eval(Const 1)) end 

= eval(reify( eval(Const 1))) 

= eval( reify(Num 1) ) 

= |val(Const 1) 

= Num 1 
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Problems of previous work 

• Double the effort: eval needs reify, print needs parse; 

• Redundancy of computation: 
eval(Appl(Abs f,e)) = eval(f(reify(eval(e)))) 

• reify must be the right inverse of eval: 

Vxe range(ev al): eval(reify(x)) = x 

eval(Abs(fn x => x)) 

= Fun( eval ° (fn x => x) ° reify ) 

= Fun( fn x => eval(reify x)) 

= Fun(fn x => x) 
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Our approach: Using an extra constructor Place to 

shortcut reify 

datatype Term = Const of int 

Appl of Term xTerm 
Abs of Term -> Term 
Place of Value 

fun eval(Const n) = Num n 

eval(Appl(f,e)) = let val Fun(g) = eval(f) in g(eval(e)) end 
eval(Abs f) = Fun( eval ° f ° Place ) 
eval(Place x) = x 

Place satisfies the desired property: eval(Place x) = x 
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Evaluating (Xx. x) 1 using the Place constructor 

eval( Appl(Abs(fn x => x),Const 1)) 

= let val Fun(g) = eval(Abs(fn x => x)) 
in g(eval(Const 1)) end 

= let val Fun(g) = Fun(fn x => eval(Place(x))) 
in g(eval(Const 1)) end 

= eval(Place( eval(Const 1))) 

= eval(Place(Num 1)) 

= Num 1 
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The catamorphism over Term 


datatype Term(a) = Const of int 

Appl of Term(a) x Term(a) 

Abs of Term(a) -» Term(a) 

Place of a 

fun cataT(c,p,a) (Const n) = c n 

cataT(c,p,a) (Appl(x,y))= p( cataT(c,p,a) x, cataT(c,p,a) y ) 
cataT(c,p,a) (Abs f) = a( cataT(c,p,a) ° f ° Place ) 
cataT(c,p,a) (Place x) = x 

fun eval x = cataT( Num, fn (Fun(f),e) => f e, Fun ) x 
eval: Term(Value) -> Value 


Restrictions 

To avoid problems: 

• The Place constructor should only be used by a cata. 

Solution: make Place and cata primitives and hide 
Place from users. 

• A contravariant variable should not be analyzed. 

eval(Abs(fn x => case x of...)) -> Fun(fn x => eval(case (Place x) of...)) 
eval(Abs(fn x => cataT(...) x)) Fun(fn x => eval(cataT(.. .)(Place x)) 

Solution: use the type system to detect these cases. 
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Doing this in general: the catamorphism over any 
algebraic type T 


It is defined in terms of the functor of T. 
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Some background 


The functor for x(a) maps a function f into the function Hx(a)](f): 
f: tj-»t 2 => |[T(a)Kf): T(ti)->x(t 2 ) 

Notation: 

fxg = X(x,y). ( f x, g y ) 
id = Xx. x 


M(f) 

= f 

[x 1 xx 2 ](f) 

= |Til(f) X [X 2 ](f) 

[T(x)](f) 

= map T ([[x](f)) 

Ix](f) 

= id 


e.g., I int x a x list(a) ](f) = id x f x map(f) 

= X(x,y,z). ( x, f(y), map(f) z ) 

Properties: [xj(id) = id [x](f) ° [x](h) = [xl( f ° h ) 





Extending functors to arrow types 


Need to differentiate positive from negative types: 
+ — + 

(int —» bool) —»int 


Positive and negative instances of a type variable a: 

+ — — + 

type T(a) = (a^a)^(a^a) 

- ' + 7 


Bifunctors 

The bifunctor of a type x(a): 

f • ti—>t2 A g: t2—>tj => |[x(a)](f,g): x(ti)—»x(t 2 ) 


|al(f,g) 

f 

[x 1 xx 2 ](f,g) = 

|x 1 ](f,g)x|x 2 l(f,g) 

[Xi->X 2 ]](f,g) = 

Alt. [x 2 ](f,g) 0 h 0 [x 1 ](g,f) 

IT(x)J(f,g) = 

map T ([x](f,g)) 

|Tl(f,g) 

id 

Properties: 

|Tj(id,id) 

S id 


M(f,g) ° [x](h,k) = [x](f 0 h,k°g) 
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The catamorphism over a type T 

Type T has n data constructors Cp Ti(T)->T, 
plus one extra constructor Place T : 


cata T (ff n ) 0 Q = ° [Ti(a)]( cata T (ff n ), Place 1 ) 

cata T (f 1? ...,f n ) 0 Place 1 = id 


datatype T((3) = C of ( T((3) 
I PlaceT of 



fun cataT(f)(C g) = f(cataT(f) ° g ° (fn h => PlaceT ° h ° cataT(f))) 
| cataT(f) (PlaceT x) = x 
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Conclusion 

Simulating the right inverse of a function with an extra 
data constructor Place is useful for: 

• meta-programming; 

• circular lists and graphs. 

It might also be useful for: 

• higher-order unification; 

• fusing programs using parametricity laws. 
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Equality between terms 


fun eqT(Const(n),Const(m)) k = 
| eqT(Appl(a,b),Appl(a’,b’) k = 
| eqT(Abs(f),Abs(g)) k 

| eqT(Place n,Place m) k 

I eqT_ 


( n=m ) 

eqT(a,a’) k A eqT(b,b’) k 
eqT(f(Place k),g(Place k)) (k+1) 
( n=m ) 


false 


eqT( Abs(fn x => Appl(x,x)), Abs(fn y => Appl(y,y))) 0 

= eqT( Appl(Place 0,Place 0), Appl(Place 0,Place 0)) 1 
= eqT(Place 0,Place 0) 1 A eqT(Place 0,Place 0) 1 
= (0=0) A (0=0) 

= true 
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Program fusion 

The Y combinator for functions: 

Y(f) x = f(Y(f)) x 

Y: Va,p. ((a^|3)^a^(3)^a^(3 

The parametricity theorem for Y is the unfold-simplify-fold law: 
Vg.Mhh: g 0 (f h) = (g ° h) => g°(Yf) = Y^ 

To find a solution for (|): 

1) choose Place such that g ° Place = id, 

2) choose k such that h = Place ° k, then: 

g ° (Y f) = Y(Xk. g ° (f (Place ° k))) 
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